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A model for the dissociation of the J/tj} is proposed, where chiral symmetry is properly imple- 
mented. Abnormal parity interactions and mesonic form factors naturally arise from the underlying 
quark sub-structure. Analytic confinement of the light quarks is obtained through an appropriate 
choice of quark interaction kernels. Dissociation cross sections of the J/t/j by either a vr or a p meson 
are then evaluated and discussed. 
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^ ; I. INTRODUCTION 

b : . . . 

, ^ , Lattice simulations of quantum chromodynamics (QCD) predict a transition from hadronic matter to a plasma 
^> ■ of quarks and gluons at a critical temperature of Tc — 175 ± 10 MeV [U, To reproduce such a condition in a 
terrestrial environment, heavy ions are collided at relativistic energies. A typical space-time evolution involves the 
collision system going through various phases including possibly the elusive quark-gluon plasma (QGP). To find an 
unambiguous signature of this new state of matter proves challenging as other stages of the fireball expansion can 
I— I. also make contributions, which will thus constitute background. 

' One popular probe of the quark-gluon plasma is the charmonium yield modification first suggested by Matsui 
^p^, and Satz in a seminal paper Q. In their original scenario, the charmonia produced in the earliest stage of the 
Q_i" collision is expected to be suppressed by the QGP due to color screening But late stage hadronic dissociation 
D could generally also occur, thus making the sources of the suppression difficult to disentangle: understanding the 
/ ■ charmonium dissociation within a hadronic gas becomes essential. However, very little is known experimentally about 
^ these hadronic dissociation processes and one has to rely on theoretical calculations, 
psj . Most studies focus on dissociation channels of the charmonia by pions as it is the most abundant particle in the 
^ produced hadronic gas. Moreover, since, in a thermal gas at realtistic temperature, the pions have just enough energy 
, to dissociate the charmonium ground-state, the J/V'j the cross section near threshold is of particular interest. Various 
' approaches can be used including non-relativistic potential models 11, B, B, Bl , QCD sum rules [ol. [Tol [ill. [l^. [isl [isl. 
ITq , and constituent-quark based formalisms HtI . fisl . flQl . [20l [2ll [2^. l23j|. Alternatively, phenomenological Lagrangians 
can also be employed 24, 25, 26, 27, 28, 29j. However, the implementation of chiral symmetry has not uniformly 
been done in all Lagrangian models [3Q] . This could then have an important phenomenological consequence, as it is 
expected to soften the cross section near the production threshold. 

In Ref . (31I . [3^ , the effect of chiral symmetry on the dissociation cross sections was investigated within a chiral 
symmetric phenomenological Lagrangian approach. It was shown that for a certain class of interactions the dissociation 
transition amplitudes should vanish as the pion momentum goes to zero (soft-pion theorem). Although, a reduction 
was observed near threshold, the effect of introducing the so-called abnormal parity interactions was at least as 
important as that of the implementation of chiral symmetry. This observation could be traced back to the fact that 
^ these abnormal parity interactions circumvent the soft-pion theorem. The results of Ref. (3^ . in particular the overall 
magnitude of the calculate cross sections, depended heavily on the parameterization of the ad-hoc form factors. It 
is the purpose of this article to address this issue by proposing a chiral-symmetric model of the J/i/j-dissociation 
whereby the form factors naturally arise from the underlying quark structure. This constituent-quark framework 
builds on models presented in Refs. p7l[l8l.lT ^ 2 ^ 21^2^[23|. However, strickly speaking, the model presented here 
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is an extension to the charm sector of Refs. (ssT 35j which itself is a generalisation of the Nambu Jona-Lasinio 
(NJL) model [s^, [s^l where the four-point interaction kernels are non-local and chiral symmetry is implemented at the 
quark level. Within this formalism, the mesonic form factors can then be calculated, and by choosing the appropriate 
light-quark kernels, it is possible to push to higher energies, or even remove, the unphysical production threshold 
of the p meson into the continuum (i.e, into a. q — q pair) which at zero temperature and density should not occur. 
This then also permits the calculation of the p-induced J /ip dissociation which is used to assess the effect of light 
resonances on the overall dissociation strength. 

This article is organized as follows: after introducing the quark interaction kernels and considering their properties, 
quark propagators are discussed with an emphasis on ways of generating quark confinement. Meson bound states are 
then found along with meson-quark vertex functions. Three- and four-point meson interactions mediated through 
quark loops are written down. Finally, after fixing the various parameters using a combination of lattice results and 
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empirical information, the behaviors of meson propagators, vertices, and cross sections are examined. The analytic 
continuation prescriptions are detailed in Appendix \^ In Appendix [Bj the isovector axial Ward identity is explicitly 
checked, while in Appendix [C] decay processes used to fix the parameters are evaluated. Finally, the amplitudes of 
the various dissociation channels studied here can be found in Appendix [Dl 



II. QUARK INTERACTION KERNELS 



In this model, chiral symmetry is implemented at the quark level. As in Ref. |32l |. the mesonic content includes 
the TT, p, J/4', D, D* , and the chiral partners of the open charmed mesons, namely Dq and Di mesons, in order 
to describe the dissociation processes of interest. Besides these, the a and ai mesons will also naturally appear in 
such an approach due to chiral symmetry. However, as for the NJL model, only global color invariance is introduced 
allowing to account for the number of quark colors, Nc, in QCD. Moreover, only color singlets are written down. 

With this in mind, the minimal action is 



S= J dx {q{x) {i^ - ml) q{x) + Q{x) {i^ - m'^) Qix)} + S„,t 



(1) 



where q and Q are the fermion fields for the light and heavy quarks, respectively, and the interactions are decomposed 
into 



with 



Y]_ / dxk 

.k=l 



Q _ C99 I cQQ 1 cQQ 
Omt — -t- + 



Kilil{x^,X2,x:,, x^Wf, {x2Wf, {X3)rf, {X4). 



The kernels, K, can be further decomposed into 



(2) 



(3) 



(4) 



where the Hi parametrise the strengths and the profiles of the non-local interactions and /; labels the flavor (either 
q or Q), T = 7or^7o. The Dirac and flavor structures for the present model are then 



{S,P,V.A} = il'«75T ,7mT ,7^757- }: 



^ V — 7m- 



Chiral symmetry then imposes that i/|''^'J'51 



H^p^^'^^^ and i/{.99^9'3> 



(5) 



j^Uii-iQ} xhis choice can be explicitly 



checked by using the global transformation for the light quark field [32l | and remembering that the heavy quark field 
is invariant under this symmetry. 

The most general form of Hs is constrained by translational invariance. To make this property explicit a change of 
variables as in Ref. [38] is made, namely 



a: = i (-a;i - X2+X3 + X4^) , 
X = X2~ xi , 

With these and suppressing all indices, the kernels become 



X' ^2 i^l + X2+ X3+ Xi) , 
x' = X4 — X3. 



(6) 



H{'f%xuX2,X3,X4) = n / '^Pfce-^^."^ P^i//^^^(pi,p2,P3,P4 



(7) 
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where the momenta are taken to be in-going. Translation invariance then amounts to requiring that under an arbitrary 
shift by a four vector a, the kernels respect 

h(^-^^{xi + a,X2+ a,X3 + a,X4 + a) = h{^^'^ {xi,X2,Xz,X4), (8) 

or more specifically that the Hi do not depend on X' . This then restricts their form in momentum-space to 

Ht^'iPi,P2,P3,P4) = i2n)U^'\P')Hl'^^p\p,P) (9) 

where we have defined 

P' = {Pi +P2 +P3+P4) 

p' = 1{P3-Pa)- (10) 
The original NJL model is given by 

h!'^'{pi,P2,P3,P4) = {27r)U^'\P') (11) 

where the G{'^^'^ are the interaction strengths, which have a dimension of inversed energy squared. Here, we will 
consider a fully-separable interaction in momentum-space inspired by an instanton-based approach [33 |. i.e., 

Ht^^{p4,P2,P3,P4) - l{27rrG{'^^S{P')ffMffAP2)ffMffM (12) 

where f{pi) are the quark form factors modeling the non-locality of the interactions. These are normalised to one 
at zero impulse and their specific forms will be chosen, in the light sector, to provide confinement. Moreover, they 
act as UV regulators for the loop integrals removing the need for a UV cutoff as in the original NJL model. The 
Gi constants scale like 1/Nc, which can be inferred by considering the simplest four-point interaction in QCD : the 
one-gluon exchange interaction between four quarks. The 1/Nc scaling of the four-point vertex will be used in what 
follows to determine an approximation scheme consistent with chiral symmetry. This particular choice of kernels also 
greatly simplifies the search for mesonic bound states and facilitates numerical integration. 



P ^ ^{Pi +P2-P3-PA) 

P^ ^{Pl~P2) 



III. QUARK PROPAGATORS 
A. Light quark sector 

The general solution of the Schwinger-Dyson equation (SDE) in momentum-space is given by (40| : 

S,{p) = Z,(p)4^^ (13) 
p^ - m^{p) 

where ruq and Zq are the momentum-dependent mass function and wavefunction renormalization respectively. 

Working in the mean field approximation or equivalently at leading order in l/Nc, the light quark propagator for 
the action of Eq. ([3]) reduces to i39,] 

S,ip) = , \ , (14) 

where the dynamical mass is given by 



mq{p) =ml + iGsfqip) J j^f,ik)TT{Sq{k)} . (15) 

We note that only the scalar channel gives a non-zero contribution and the dynamical mass scales like as in QCD 
[39| . The above gap equation then admits the solution [33] 

mqip) = ml + imqiO)-ml)f^{p) (16) 
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where mg(0) is the dynamical mass at zero momentum. The quark propagator can also be directly linked to the quark 
condensate in the chiral limit through the expression 

(99>o = -^^TrlS.ik)]. (17) 

From the quark propagator it can be readily inferred that there will be no poles on the real axis in two cases: either 
the poles are complex [19| or there are no poles at all in the complex plane Here, we will consider the second 
case. However, to illustrate how this property is manifest we need to analytically continue the quark propagator to 
Euclidean space, i.e., the denominator becomes p\ + 'm?{p). 

We then follow Ref. |43| where the inverse of the quark propagator denominator, in the chiral limit, is parametrized 

as 

1 1 - e^^P^ 

(18) 



p2 -(- 777,2 (jy^ p2 

Alternatives also exist such as the one found in Ref. [i^. Various limits can then be considered. For large positive p^, 
this quantity behaves as expected perturbatively, i.e., 1/p^. This is not the case for large negative p^ where it diverges. 
This is not a major problem provided the mass function is probed only for small negative p^, i.e., small time-like 
separation in Minkowski space. As the infrared limit is approach, i.e., p^ —^ 0, the inverse quark propagator becomes 
constant and equal to /i. The cutoff parameter is then seen to be equal to fi — l/m^(0). This is sufficient to show 
that no poles exist. Mathematically, this function property is called entire, and this realization of quark confinement 
is deemed analytic. Reinstating the current quark mass yields 



mq{p) 



(19) 



where the principal branch has been chosen.^ In what follows, we will refer to the model using this functional form 
for the dynamical mass as Model I (MI) . This choice will complicate the evaluation of n-point functions as branch 
cuts will generally appear. Since the zero-momentum light-quark mass should be of the order of a few hundred MeVs 
0, H^li the parameter ^ will be greater than one. Therefore, for p^ > 0, the dynamical mass will exhibit a steep 
decrease. To control this ultraviolet behavior, we also consider an alternative dynamical mass model (Mil), which is 
constructed by substituing — > + (mo — m^)e""P in Eq. (jl9[) where toq and a are two additional parameters. 
The steep decline for > can then be overcome. 



B. Heavy quark sector 

For the heavy flavor sector, again no wavefunction renormalization is possible for this model at leading order in 
l/Nc- Furthermore, there is no heavy quark scalar interaction in the quark action. Therefore, there is no dynamical 
mass generation at the mean field level contrary to the light quark sector. This leaves the charm quark form factor 
unspecified. Here, we will chose: 

/Qb) = e-^«P' (20) 
in Euclidean space where (3q will be fixed by considering the decay of the J/ tp into dileptons (see Appendix [C| . 



IV. MESON BOUND STATES 



A. Meson-quark vertex functions and meson propagators 



Interactions amongst constituents lead to the emergence of bound states provided attractive channels exist between 
them. In QCD, these occur, at least perturbatively, in the singlet-color channel between a quark and anti-quark. At 



^ This choice, through Eo. 1161 1. leads to an action that is not linear in the current mass, contrary to the QCD action, and results from 
enforcing analytic confinement in this particular model. However, it can be shown that the GMOR relation holds to first order in rric 
for such a mass model (see Ref. [42| ) . 
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all orders in quantum field theory, bound states are found by considering the pole structure of the scattering matrix, 
S. Following Refs. 0, E^, the inhomogeneous Bethe-Salpeter equation (BSE) satisfied by the T-matrix for a given 
flavor content, i.e., the interacting part of the scattering matrix, is 

T{p,p',P) = K{p,p',P) + I dp"K{p,p",P)S{p-)T{p",p',P)S{p+) (21) 

where the Dirac indices and flavor labels have been suppressed for clarity and p± — p" ± ip. 

To solve the BSE, we have to specify, besides the quark propagators, the scattering kernels. In general, since they 
comprise of an inflnite sum of diagrams, a truncation scheme must be implemented. If the Hartree propagator for 
light quark is used and we require chiral symmetry to be maintained, then the scattering kernels involving light quarks 
are uniquely determined: they cannot contain exchange terms and their functional forms are fixed to that of Eq. (|12|) . 
This is known in low-energy nuclear physics as the random-phase approximation (RPA) and is the counterpart of 
the mean field hypothesis for solving the SDE. Moreover, without any surprises, their iVc-scaling is that of the Gi 
couplings. The kernels being fully separable, the BSE admits the solution 

T^^^^ ip,p',p) = ffApi)fhiP2)ffAP3)ffMff^^HP) (22) 

where momentum conservation is implicit and T^'^^^{P) are the T-matrix with the f{p) factored out. 

Next following Ref. [46i|, the BSE is decomposed into independent Dirac channels. This can be achieved by using 
the projection operators 

^t^iu — 9^1^ P^Pu-j Lf.Lv — PfjiPf (23) 
where P^ = The scattering kernels can then be re-arranged for the pseudo-scalar and scalar channels as 

)J'/i/2 /^/i/2 / p/1/2 ^ p/1/2 \ 1/^/1/2 l'nfif2;L„~pfif2;L\ (OA') 

^{S,P} - ^{S,P} 1,^ {S,P} ^ {S,P} J + '^{V,A} \^ {V,A} ^ {V,A} J ^ > 

and 



T>f\fl _ r^S^fl ( -p/1/2 ;T „ p/1/2 ;T\ lo^\ 
^{V,A} — ^{V,A} Y {V^A} ^ ^ {V,A} ) y'^^l 



for the vector and axial ones where 



^ {V,A} - + ^v)^{V,A} - ^ {V,A} l^"-* 

Note also that the left- and right-hand parts of a given fiavor and Dirac matrix product is associated with the in- 
and out-going quark states, respectively. Thus the kernels should be read from left to right. Similarly, the T-matrix 
can be decomposed into products of Lorentz covariant tensors giving 

7^/1/2 _ 71/7-/1/2 /-p/1/2 ^ p/1/2 A I 7i/r/l/2 / pi/1/2 ^p/l/2;i 

{s,p} ~ ^^^{ss.pp} y- {s,p} ^ ^ {s.p} ) ^^^{sv.PA} y- {s.p}^ ^ {v,A} 

I Af/1/2 (yfih-.L yfiS-i \ _i_ ]\,ffih\L (yfif2;L j^fif2;L\ /r,yN 

^ ^"{ys,AP} {v,A} ^ ^ {s,p} J + ^"{yy.AA} {v,a} ^ ^ {v,a} ) ^ y^^' J 



and 



71/1/2 _ i\/rf'^h\T 

{V,A} — ^^^{VV,AA} 



where the M components are functions of P^ . It is clear that, in general, mixing occurs between channels, i.e., the 
longitudinal component of the vector and axial channels contribute to the scalar and pseudo-scalar ones respectively, 
and that the T-matrix is then block-diagonal. From the BSE, it is then seen that for a given flavor and Dirac 
channels, a given M matrix satisfies the equation 

M = G[l + JM] = ^^ (29) 

where all labels have been suppressed. The fermion-loop matrix is then given by 



ijl^f^{P') = - / dkfl{k,)fl{k+)Tv rf ^^%(ft_)rf ^^%(fe+) (30) 
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where the minus sign on LHS is due the normal ordering of fermion fields, k± = ± -j, and the matrix indices are i 
and j.^ This equation shows that a pole will develop in the T-matrix when 

A=\l-GJ{mli)\^0 (31) 

where mj\/ is the mass of the meson. On one hand, near this point, the T-matrix admits the solution 

Tip^p'^p^^iMEin^i^^iip^in (32) 

P rrij^j 

where the wavcfunction for meson M is defined as 

xm{p', p) = (97w/. \M{P)) = 9Mff, {Pi)fh iP2) f 1 - —r) r{/^ (33) 

with p' = 5 {ps ^ P4)- The coupling and mixing parameters are found by solving Eq. (|29[) using Ea. (|32p yielding for 
a meson of spin s: gm — {—lyM'^^/dA/dP'^ and = M^^/M'^^. Note that for channels where there is no mixing 
ajvf = 0. On the other hand, for far from the on-shell condition, the T-matrix can be written as 

Tm = igM^M iguTiPM {P^ ) (34) 
where T\j = (l; ® ^{^/^ ^^'^ with the meson propagator given by 

VMiP) ^ ~^Mm{P) (35) 
9m 

where the tensorial structure is suppressed. 



B. Meson interactions 



Having discussed the couplings between quarks and mesons, the interactions amongst the mesons are now examined. 
As for the meson self-energies, only interactions mediated by quarks will be studied; mesonic fluctuations being sub- 
leading in a l/Nc expansion. 




P ^ 



P k -P 



k+ P 




FIG. 1: Momentum conventions for the three-point vertex functions. 



To evaluate the dissociation cross sections, the three- and four-point interactions have to be written down. The 
former are further divided into two, namely interactions between one light meson, either the pion or the p meson and 
two open charmed mesons; or interactions between the J/ip and two open charmed mesons. The momentum flows for 
these two cases are depicted in Fig. [TJ Here, in general, one of the three mesons will be off-shell and the kinematical 



^ The analytic continuation for these two-loop integrals is explained in Appendix 1X1 
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variable t is then the associated momentum transfer. The expressions for the meson form factors are 

dfcTr 



dkTr 



iXLiSq {k-)igHXH.,iSQ (fc+ - PffJ iXH.iSq (fc+) 
iXLiSq ixH2-iSQ (fc+ - PgJ igHit' Hi'iSq (fc+) 



(36) 
(37) 



and 



dkTi 



dkTi 



iX^iSq (fc - Ph^) igHj^H3^Sq{k)ixH:,iSQ {k + Ph^) 
iXi'-iSg {k - PgJ ixH3iSq{k)igHj^H2iSQ {k + Ph^) 



(38) 



(39) 



where the momentum arguments of the wavefunctions have been suppressed for clarity, L and H label the light and 
open charmed mesons respectively, the minus sign is due to the fermion-loop, and the round parentheses indicate 
which meson is off-shell. The extra label i notes that the form factor can be a two-component vector due to mixing. 




k + P 



FIG. 2: Momentum convention for the four-point vertex functions. 



Similarly, the general expression for the four-point interaction, illustrated in Fig. [21 is given by 



i XLiSq{k_)ixH3iSQ{k^ + PH3)iXi^ 
i Sq {k+ - PffJ ixaJSq (fc+)] 



(40) 



where again the minus sign is due to the fermion-loop. The analytic continuation of the loop integrals is discussed in 
Appendix [XI 



RESULTS 



A. Parameter fixing 



The isospin-averaged masses of the pion, p, D, D* , and J/tp mesons used in what follows are 0.138 GeV, 0.770 GeV, 
1.868 GeV, 2.009 GeV, and 3.096 GeV respectively. The first model (referred as MI) is that of Eg. p9)) and has four 
parameters, namely m^, /3q, /it, and m^. Three observables, i.e, the pion decay constant [Eq. (|C2ip ]. the dielectron 
decay width of the J/ip meson [Eq. (jC6|) ]. and the the decay width of the D* into D-tt final state [Eq. (|C3[) ]. can be 
then used to partly constrain the possible parameter values. Specifically, ml is determined in order to reproduce the 
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While, for a fixed m^, the 
This leaves the 



pion decay constant value since its expression [Eq. (|C21|1 ] has a strong m^-dependence 
range parameter, /3q, is constrained by the experimental dielectron decay width of the J/V' [Ea. (|C6l 
open charmed decay constant [Eq. (jC3j) ] to fix both and fj,. To restrict the possible solutions, we also required that 
the light-quark dynamical mass at zero-momentum be of the order of a few hundred MeVs as it is seen from lattice 
results 38] and as used in other phenomenological studies 0, H^l • We further require that obey the constraint 
mj/^ < 2m^ since the model is non-confining in the heavy sector. With these additional constraints, one possible set 



of values is = 1.59 GeV, /3q = 0.06 GeV" 



12MeV, and fi = 8 GeV"''. These then yield = 94.78 MeV, 



-+e- 



J/4!'^e++e- 



5.44 keV, and Qt^dd- = 18.62 which are to be compared to the experimental values of — 93 MeV, 
5.5 ± 0.14 ± 0.02 keV and g^^oo, = 17.9 ± 0.3 ± 1.9 [48]. The resulting light-quark mass at zero 



momentum is then mq(G) — 0.354 GeV, while the quark condensate is —(237)"^ MeV as calculated using Eq. (fT7)) . 
Both these values are slightly higher than those found in Ref . [s^] . 

Introducing the substitution discussed earlier in Section llll Al into the Eq. the number of parameters in the light 
sector increases by two (model Mil). Additional information is thus needed. An alternative could be to calculate other 
observables such as the p decay into two pions in order to fix the extra parameters. Rather, here a parametrisation 
of lattice data is used.'^ Specifically, we will utilize the one proposed in Ref. [38j] , namely 



ruq (p) 



A2 



A2 



(41) 



where am = 0.343 GeV and A„j — 0.767 GeV. Fitting our model [Eq. (fT9|) ] to this parametrisation for e [0, 1] gives 
Too = 0.227 GeV, a = 1.096 GeV"^ and ^ = 11.786 GeV~^. The light current mass is then taken to be = 5.5 MeV 
yielding /^r = 92.11 MeV. The zero-momentum dynamical mass is then 0.331 GeV, while the quark condensate is now 
-(239)3 MeV^ With the same values for the heavy-sector parameters, the coupling constant g-nOD- is 18.35 which is 
a prediction of the model and well within the experimental tolerance [481 ]. 




0.5 

p**2( GeV**2 ) 



FIG. 3: (Color online) Light -quark dynamical mass models with u = p^. The solid and dashed lines are the MI and Mil models 
respectively. 

Fig. [3] shows a comparison between the two light-quark mass models discussed. As seen, for < 0, the two 
models are very similar. The main difference is for p^ > where for MI the dynamical mass drops very quickly down 



^ The quark mass is not a gauge-invariant quantity and thus dependent on the gauge chosen to carry out the simulation. The purpose 
here is to capture some flavor of the QCD dynamics. 
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to its asymptotic value of m^. Both models have similar zero- momentum mass, quark condensate, and calculated 
observables. The only noted difference is in the light-quark current mass for which there is more than a factor two 
difference between models. Although the current mass for Mil is within the range given by PDG [i^, i.e., 1.3 to 5 
MeV for the u-quark and of 3 to 7 MeV for the d-quark, one could wonder if introducing two extra parameters, toq 
and a, in order to reduce the current quark mass value is justified. At this point, we will retain the two models to 
ascertain if any other differences occur for vertices and cross sections. 

With the parameters for MI and Mil, the four-quark couplings, the meson-quark couplings, and the mixing coeffi- 
cients can be evaluated. Their values are listed in Table HI Since the Dq and Di are the chiral partners of the D and 
D* respectively, the masses are not independent and have to be calculated by finding the zeros of the respective me- 
son propagator denominators. Doing so yields m^). = 2.064 GeV and = 2.249 GeV, and, m^)* — 2.045 GeV and 
rriD^ = 2.231 GeV for MI and Mil, respectively. These are to be compared to the experimental masses of to^i* = 2.40 
GeV and mn^ = 2.43 GeV. It is quite clear that neither model is capable of reproducing the absolute masses and the 
mass difference, i.e., Arriexp = 0.03 GeV. This problem then implies that the interaction kernels or, in non-rclativistic 
terms the potentials, require further modeling. This is left for a future study. 



Gm 


MI 


Mil 


9M 


MI 


Mil 




-1.145 


-1.145 


9J/ip 


1.717 


1.717 


Gd' 


-6.690 


-5.257 


go* 


2.025 


1.842 


Gdi 


-6.690 


-5.257 




1.955 


1.772 


Gd 


17.661 


11.977 




4.667(0.301) 


4.166(0.315) 




17.661 


11.977 




3.654(0.166) 


3.828(0.205) 


Gp 


-7.070 


-6.147 


gp 


1.336 


1.219 




52.562 


31.052 




3.768(0.0220) 


3.615(0.0233) 



TABLE I: Quark-quark couplings, meson-quark couplings, and mixing coefficients. 



B. Meson propagators and vertices 

Having fixed the parameters for the two models, the meson propagators and vertices can now be examined. For 
the meson propagators, I?m, as can be seen from Eq. (j29[) . the asymptotic behavior is controlled by the the two-point 
functions, Jm, as the quark-quark coupling, Gm, is independent of momentum.* As t — > cxd, Jm — > oo and Dm ^ 0, 
while for i — > — oo, Jm and T>m ^ — ^r^. Near the meson pole, the scalar part of the corresponding propagator 

[Eq. ((32)) ] is expected to behave as 

where s is the meson spin. Note that this form is used for the phenomenological Lagrangian studies of Refs [l^, [H, 

Fig. m ^ illustrates the behaviors of the D* and p propagators for MI and Mil compared to that of Eq. ([1^ . It is 
quite clear that the phenomenological propagators are comparable to the NJL ones only near the poles, and that the 
differences between the propagators of MI and Mil are not significant. Although not shown on the figure, this latter 
observation holds true for t > m\j. Note that for t ^ —oo the NJL-propagators exhibit the expected asymptotic 
behavior, i.e., it is non-zero as for the Lagrangian one, but rather — t"-. 

9m 

Fig. [5] shows examples of three-point meson vertices. The curves labeled Lagrangian are the ad-hoc form fac- 
tors multiplied by the relevant meson couplings used in Ref. [s^]. Let us first consider the differences between the 
Lagrangian and the NJL approaches. We note that a relative agreement only exists for the coefficient of ^Vd(D') prO" 
portional to pion momentum, while the overall magnitude of the coefficient proportional to the p-meson momentum 



For simplicity, mixing is ignored here. Adding it does not alter the conclusions. 
^ A logarithmic scale was used in order to permit the discrimination between the three curves. At the pole, the divergence should be 
infinite and the appearance of the contrary is just an artifact of the finite number of points in the numerical evaluation of the propagator. 
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FIG. 4: (Color online) Scalar parts of the transverse p and D* meson propagators near their respective pole. 



of F^Q^j-^^ is suppressed and the energy dependence of the coefficient proportional to the J/i/' momentum of F'j/^Di^Q-j 
is quite different. 

Turning to the comparison between MI and Mil, we note that the vertices are quite similar both in their overall 
magnitude and their energy-dependence with the exception of the ploted coefficient of F^Q(^ jj,y For large space-like 
separation, large fluctuations appear for MIL This can be linked to the difficulty of carrying out the principal value 
integral in Eq. ([39|) due to the large oscillations in the heavy quark form factor in the loop integral. As seen in Fig. [5l 
increasing the number of numerical evaluations reduces the fluctuations. Although, not apparent in the figure, this 
problem exists for all vertices with a light meson. 

Finally, the four-point coefficients of the form factors F^j/^dd* ^'^d Fpj^^Qjj proportional to the metric tensor 
are plotted in Fig. [5] for three center-of-mass energies. We note that the form factor used in the Lagrangian approach 
of Ref. [Ill is very suppressed and quite flat compared to the NJL equivalent. Moreover, the differences between the 
form factors for MI and Mil are slightly more pronounced than observed for the three-point functions at least in 
terms of magnitude. 



C. Cross sections 



The transition amplitudes can be found in Appendix [Dl The cross sections are plotted in Figs [7] and [51 ^'^ We first 
note that the results for both MI and Mil are very similar for small ^/s and differ only slightly in magnitude for larger 
values, and therefore the introduction of two additional parameters in Mil is probably not justified. Comparison 
with the results of the phenomenological Lagrangian study of Ref. ^] shows significant differences with MI from 



Because of mixing between open charmed mesons subtle cancellations occur between sub-amplitudes. For the process tt -|- J/i/) — > 
(D -|- D*) + (D + -D*), for example, this happens between the D- and Di-exchange channels, and between the D*— and Dg-exchange 
channels. The requirement that they cancel can be traced back to the tensorial decomposition used to split the transverse and longitudinal 
parts of the vector and axial- vector Dirac structures Eq. H26| l. This splitting then induces two factors where q is the momentum 

flowing through the propagator. One is absorbed in the transverse projector of the vector particle propagator, while the other is 
further split between the two vertex functions sandwiching the spin— propagator. As — > 0, divergences appear. Analytically, when 
all the sub-amplitudes are summed, they cancel; the splitting being artificial. However, these cancellations amount to subtracting 
large numbers. This lead to a numerical integration problem. Indeed, for certain ^/s, the quadrature method employed can require 
evaluations at points close to q^ = 0. To deal with this problem, we force the cancellations within a small radius centered around q^ = 0. 
In some sense, this is an estimation of the numerical precision associated with the evaluations of the vertices and the propagators. The 
more precise the evaluations are, the better the cancellation is. A radius value of 0.05 GeV^ is used here. Relics of the incomplete 
cancellations usually still remain in the form of small bumps in the data as seen in Fig. [7] around y/s = 4.4 GeV and ,/s = 4.8 GeV for 
both -K + J/i/i ^ (D + D*) -I- (D -I- D*) and tt -|- J/ip ^ D* + D* {tt + J/ii D + D does not have divergences). 

In order to assess the effect of the finite current mass, the cross section for ■K + J/'ip^D + D was re-evaluated with a zero current mass. 
It was found that the deviation is small and no greater than 3.2% at 5 GeV and less than 1.5% near threshold. 
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FIG. 5: (Color online) Examples of three-point vertices. 



the onset. We note that the maxima of the pion dissociation cross section for tt + J/ip ^ {D + D*) + {D + D*) 
and TT + J/ip ^ D* + D* are smaller by about 50% than those found in the previous non-local NJL study ^] and 
in the potential model approach * Since in the potential model, the spin-orbit interaction is not modeled the 
TT + J/ip ^ D + D process is not evaluated Q , and a comparison can be made only with Ref . (23. Contrary to the 
two other pion- induced dissociations the maximum for this process is comparable to that of Ref. [22| . 

A consequence of the mixing between channels is that, it is impossible to consider normal parity sub- amplitudes on 
their own in order to asset the effect of chiral symmetry since doing so would entail divergences appearing. However, 
this does not prevent us from considering the effect of the abnormal sub-amplitude outt + J/^^ {D + D*) + {D + D*) 
process. In Fig. [9l we observe a reduction of the cross section near threshold when the abnormal contribution [A^26 
of Appendix [P] is removed. Specifically, the maximum {^/s w 3.92 GeV) is seen to decrease by 14%, which is far less 
than what was found in phcnomcnological Lagrangian approach of Ref. [s^l . 

For the dissociation by the p meson, we estimate a maximal cross section for p + J/tp {D + D*) + {D + D*) 
comparable to what is found in the potential model while we find a maximum for p + J/tp — > D* + D* which is 
an order of magnitude larger. Finally, the trends of the energy behaviors of all p-induced dissociations are similar to 
those found in Ref. 



In Refs I22II and specific charged channels arc plotted, while here we present isospin averaged cross sections. In order to make contact 
with these studies, the 7r(p) -|- J/-)/) (D + D*) + (D + D*) cross section has to be divided by two. 
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FIG. 7: (Color online) cross sections for the J/-!/)-dissociation by a pion. Right panel is a comparison between MI and Mil, 
while the left panel is one between MI and the results from the phenomenological Lagrangian study of Ref. (32jl for A = 1 GeV. 



VI. CONCLUSION AND OUTLOOK 

In Ref. [S^l , the absolute values of the strength of tt- and p-induced dissociations depended on the choice of form 
factors and the techniques used to fix their absolute normalizations, which put into question the robustness of the 
model. To address this problem, a non-local NJL similar to that of Refs [ij, H^l was presented. The vertex form 
factors were then calculated from the underlying quark structure, thus reducing some of the uncertainties found in 
the previous study of Ref. [Ill . We further utilized the fact that in the non-local version, quark form factors can be 
chosen in order to confined analytically the light-quark propagator, i.e., the light quark propagator then has no poles. 
Doing so permitted us to calculate the dissociation cross sections by a p meson. However, it was impossible to asset 
the effect of chiral symmetry since mixing between channels prevented the removal of the sub-amplitudes where chiral 
partners are exchanged. This problem did not affect the abnormal parity term for tt -I- J/-0 ^ iP ^ D*) + {D + D*). 
Turning it off then lead to a reduction near threshold as expected. But this decrease was far less substantial than the 
one observed in the phenomenological Lagrangian approach of Ref. [33 | . 

Further work could include calculating the dissociation of higher charmonia as it is relevant within the context of 
sequential charmonia melting in the QGP [soj . This would entail improving the heavy-quark four-point interaction 
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FIG. 9: (Color online) Effect of abnormal parity content on the n + J/ip ^ {D + D") + [D* + D) cross section. The dashed 
and solid lines correspond to cross sections without and with the abnormal parity sub-amplitude. 



kernels in such a way that they could sustain higher resonances. Similarly, other light mesons and, potentially, open 
strange mesons should be considered. Finally, calculating semi-leptonic decay constants and other observables, such 
as the electromagnetic pion form factor, could lead to an improvement in the modeling, which would then increase 
the confidence in the overall magnitude of J/T/i-dissociation cross sections. 
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APPENDIX A: ANALYTIC CONTINUATION 

In order to calculate the various loop integrals, an analytic continuation to Euclidean space is needed. In pertur- 
bation theory, the usual prescription is the Wick rotation, i.e., for the loop momentum: k'^ —^ ik^, which ensures that 
the Feynman boundary condition is encoded and that only the relevant poles are picked up [52,]. For a model with 
analytic confinement, a problem arises. This can be made explicit by examining an example. 

Consider the loop integral of Eq. ([50]) for the self-energy of a light meson after a Wick rotation 

tJ,^{P'^) ^-ij dkE fl (kE ~ flP) fl {{kE + (1 - ri)P) 

X Tr [r^Sg {kE - vP) i'jSq {{kE + (1 - v)P)] (Al) 

where an arbitrary momentum shift of the loop momentum parametrised by 7] has been done and P = (^—iP^, P^ with 

P^ — —P^. Note that the squared light-quark form factor has a square-root dependency since it is defined through 
Eq. (|16p in the Hartree approximation. Furthermore, in Euclidean space, the arguments of the form factors and the 
quark propagators are in general complex. The conjunction of these two elements can then lead to the appearance of 
branch cuts making the integrand ill-behaved. 

To make this explicit, the chiral limit of the loop integral is taken, namely 



JiP%j=- I dkE J^- ^ g{k%,kE-P,P') 

V -2— rj-^l — ri 



(A2) 



where the arguments are defined as zs — {kE + SP)^ and the function g {k'^, kE ■ P, P^) is the product of the Dirac 
and flavor traces. An Euclidean coordinate transformation is then performed to go from the Cartesian coordinate 
system to the spherical coordinate one. This is done through (53 | 



kE = {k4,k) — > {ky, k\J \ -ip- \/ \ — cos (j),k\/l — y^\/\ — sin 4>,k\/l — y'^x). (A3) 

Furthermore, to reduce the number of non-trivial integrals to carry out the fermion-loop is evaluated in the rest-frame 
of the meson. The complex arguments become 

z_r, ^ k^ + 2irikyP° -r]^P^, zi-^ ^ k^ - 2i{l - jj)kyP° - {1 - rj)'^ P^ (A4) 

with P — {mM, 0). 

For the case of P^ ~ m^, both the real and imaginary parts of the integrand with the trace function g factored out^ 
are plotted as a function of k in Fig. [10] for y = 0.9 and two values of r]. We note that for 77 = 0.1 the imaginary part 
has several discontinuities which occur when the real part is zero. However, for 77 = 0.45, both the real and imaginary 
parts are smooth functions of k. This latter statement turns out to be true for all y in [—1, 1]. Moreover, such a 
behavior can be identified for a range of 77 values. It is then possible to evaluate the integral for several rys within this 
range and check that the results are equal as expected from translational invariance. 

We note that by choosing r] — \ the arguments of the square-root becomes real and positive for all value of k and 
y. In order words, the two squared light-quark form factors are complex conjugate of each other. It is important to 
remark that this is the case only because the evaluation of the fermion-loop is carried in the rest-frame of the meson. 
Thus by doing the appropriate shift of momentum, the evaluation of the integral in the light-meson's rest-frame 
is numerically tractable. Moreover, this technique can be applied straight-forwardly to loop integrals of three- and 
four-point correlations. 



An overall minus sign and the constant numerical factor have also been omitted. 
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FIG. 10: Real and imaginary parts of the integrand with the trace function omitted. The solid and dash curves are the real 
and imaginary parts respectively. The left panel is for 77 = 0.1 while the right panel is for rj = 0.45. 



However, the above procedure cannot be used when evaluating correlation functions with no external light-meson. 
The simplest example is for the open charmed self-energies where there is only one squared light-quark form factor 
rather than two. The requirement that there are no branch cuts can then be implemented by shifting the loop 
momentum in such a way that the integral reads 



(A5) 



The argument of the squared light-quark form factor is then always positive and real. However, doing so does leave 
the possibility that for a certain the Euclidean heavy quark propagator can go on its mass-shell for some of the 
fc^j-integral points. The required continuation prescription can then be found by going back to the Feynman boundary 
condition. The heavy quark denominator in the loop integral can be written in Euclidean Cartesian coordinates as 



(kE +P)+ 



k/\ 



2ik4^mM — Tom + 



Q 



(A6) 



where niM is the meson mass. When 



— and P^ — 0, the poles are at ±imQ and the Feynman prescription 



dictates that only the positive pole residue should contribute to the line integral. Keeping the meson mass to zero. 



but increasing 



, we note that the poles move away from the origin. Thus the contour depicted in upper-left panel 



of Fig.[TT]is equivalent to the Wick continuation. Reinstating a meson mass leads to two cases. The first one is when 



mj\/ < mg . For a null loop three-momentum, the poles are on each side of the real fc4-axis and as 



increases they 



move away from each other. In this case, there is again no need to alter the contour. The second case, which is more 

interesting, occurs when uih > mq. Both poles are in the upper part of the complex plane for k = 0. One of the 

pole eventually migrates to the lower half-plane as the loop three-momentum is increased. As this pole crosses the 
real A:4-axis a jump occurs in the line integral. This is due to the fact that by evaluating the line integral along the 



= should not contribute to the integral, 



fc4-axis, the residue of the pole, which in the limit where ijim = and 
is included. The solution is then to deform the contour as in the lower-right panel of Fig. [11] to exclude this pole 
between 



= and 



m 



M 



Iq where the latter point is found under the condition fc4 = 0, i.e., when the 

pole is on the real fc4-axis and is about to go into the lower half-plane. 

The above example is one of two possible scenarios generally encountered. The other one happens when P — > — P 
in Eq. I|A5|) . For tuh > mq both poles start in the lower half-plane. Thus, the contour has to be deformed now to 

include the pole required by the Feynman boundary condition for the k -interval where it is in the lower half-plane. 
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(not shown here). 



This description can be systematically extended to higher-point correlation functions. The final integration pre- 
scription in Euclidean space is then, after choosing the appropriate loop momentum flow, to evaluate the principal 
value of the line integral and add or subtract the appropriate residues. 



APPENDIX B: ISOVECTOR AXIAL WARD IDENTITY 



It is important to verify that the approximation schemes for the quark and meson propagators are consistent with 
each other and do not break chiral symmetry in the chiral limit. In our model, the divergences of local currents are not 
zero. Rather residual terms due to the non-local interactions are left. This problem of constructing a gauge-invariant 
non-local theory is well studied and we refer the interested reader to the Refs [3l,|33,|5l[ for a complete treatment. 
In particular, in Ref. [s^] the vector current is explicitly constructed and the related Ward identity is checked. Thus, 
only the isovector axial Ward identity has to be ascertained in order to ensure that chiral symmetry is valid within 
the approximation context. 
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For the isovector axial symmetry, its Ward identity is [4C 



(Bl) 



where the momentum flows are given in Fig. [T^ ^'^ Two cases are considered. The first one is when there are 







FIG. 12: Contributions to the effective isovector axial vertex. The diagram (a) is the local contribution while diagrams (b), 
(c), and (d) are the non-local contributions. 



only scalar and pseudo-scalar four-quark couplings, while the second one includes mixing due to the introduction of 
a vector and axial channels. In both cases, the term due to the current quark mass is omitted. For the first case, 
they are four different contributions to the isovector axial vertex. These are depicted in Fig. [121 The first non-local 
contribution to the divergence (diagram labeled b in Fig. I12p is due to a scalar fermion-loop which can be inferred 
from Eq. (13) of Ref. [s^l- Summing the local and first non-local terms yields 



^t^^ 5 



(B2) 



where the scalar fermion-loop is defined as 

Is{P) = Gs j dkfg{k)TT [Sg{k)] (/,(fc + P) + f,{k - P)) . (B3) 

Next, the contributions due to the pionic resonance in the absence of mixing [diagrams (c) and (d) of Fig. I12j can 
be cast as 



.56 _ . Gpfq{pi)fq{p2) ^^^^_^a 



^\-GpJpp{P) 

Inserting Eq. (|B2p into the above equation yields the final expression for 



dkfg{k+)fg{k^)Tr j5r''Sg{k+)P^rfSg{k 



(B4) 



r5fc _ _ • 



1 - GpJpp{P) 2 ^ F FF\ )\ b\ ) 



-^h{Pi)fMT.^Is{P). 



(B5) 
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From Eq. I I12I I. it is expected that for a finite current mass; P^Fg'' = 2mc.75^. 
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Summing this contribution and that of Eq. (IB2[) verifies the axial Ward identity. 

Adding the vector and axial channels leads to an additional contribution to Eq. (|B2p due to the vector insertion in 
the fermion-loop. The divergence of the resulting non-local current in momentum-space is inferred from Eq.(lO) of 
Ref. \3M and reads 



/p2 '""2 



where the vector fermion-loop is 

Iv{P) j dkTi -^Sqik) 

The contribution due to the pion intermediate state is then 

r% = -tfMfM{isiP) + 



U{k) {Uk + p)-Uk-p)). 



Iv{P)\l5 



(B6) 



(B7) 



(B8) 



which again cancels both the second term of Eq. (|B2p and Eq. (jB6|) thus verifying the axial Ward identity for this 
extension. 



APPENDIX C: DECAYS 
1- gTvD'D coupling 

For the dissociation processes studied, all three-point vertices are evaluated with one external particle off-shell. 
Moreover, for most of them the kinematics do not allow to have all three mesons on-shell. One exception is for the 
D* ^ D + TT decay process which has been measured experimentally by CLEO |48| 

We then wish to use this experimental information to constrain the parameter set. To do so, the expression of the 
decay width as a function of the on-shell three-point coupling is written down: 



y-nD-'D 



p. 



48m|,. 



(CI) 



where 



P. 



is the centre-of-mass three-momentum. The second step then consists in relating the coupling to the the 



extended NJL model. The associated meson form factor can be decomposed into 



^ D*+- 



with the coupling given by 



gnD'D=V2{Fo-Fi) 



(C2) 



(C3) 



where four-momentum conservation Pu* = P^ + P^ and orthogonality e (Pd* ) • Pd' — have been used, and the 



factor is to account for isospin ( the coupling gj^D* d being defined to be equal to gjj. 



|48l). 



2. J/tp decay into a dilepton 

We now turn to the calculation of the decay of the J/tp into a dilepton. Fig. [T3l shows the contributions to the 
effective quark-photon vertex. It is assumed that the dominant behavior will be given by diagrams at leading order 
in 1/Nc- Thus, only the transition due to constituent quark loop in the direct channel will be considered; exchange 
diagrams and mesonic fluctuations are ignored. It is important to note that this approximation is consistent with the 
Ward identity [sj ■ The transition amplitude of an on-shell J/ifi into a photon is then 

zA^^;(P) - - j dkTT[zx^{k,P)zSQ{k^)iTl,,,{k,P)zSQ{k+)] 

= -g,peQT^i j dkfQ{k+)fQ{k-)Ti-[rSQ{k-)T'^,j,j{k,P)SQ{k+)] 

= -g^eQT'^''[l^{P)+I^'^iP)] (C4) 
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where Q is the electric charge of the heavy quark, e = |e|, the minus sign in front of the integral is due to the 
fermion-loop, and the wavefunction meson label is suppressed. The coupling constant between the full current and 
the photon field, is derived from 



Cem = ~eQJ'^j^^{x)A^{x). 
And coupling the photon to a dilepton finally yields the decay width 



r 



e 



4 



\P\- 



(C5) 



(C6) 



where |Pp — ^ — mf. 

The main ingredient left to specify is thus the effective electromagnetic vertex. In the heavy quark sector, only a 
vector coupling is introduced. Thus there is only one correction term to the local vertex and we can write 



where 



and 



Il^iP) = - f d^j dkTr[j,SQ{k)]-^fQ{k + XP)fQ{k-{l 



Inserting the electromagnetic vertex into the transition amplitude gives 



^M^•'(P) 



and 



-(l-A)P). 



-g^eQn j dkfQ{k+)fQ{k^)TT [7"^Q(fc_)7^^Q(fc+)] 



(C7) 
(C8) 

(C9) 



where, for the non-local term, Eqs. ([50)1 and (pij) have been used and the scalar integrals are defined as 

4^C / dkfQ{k+)fQ{k^) 



2 {k-py 

3 



kl- 



(ClO) 

(cii) 

(C12) 



and 



tNL _ 

Irp — 



d\ / dk 



p2 



fc2 



dk' 



jQ{k + XP)fQ{k~{l~\)P). 



(C13) 



From the non-local scalar integrals, we note that the final result will dependent on the interpolation path. This is 
due to the fact that the coupling between the photon and the J/^ is transverse and, thus, not constrained by the 
underlying current conservation. 



20 



3. Pion decay constant 

The leptonic decay of the pion into a muon and a muonic anti-ncutrino can be studied by considering the couphng 
of the pion-field to the isovector axial current [SSj . Formally, the coupling is inferred from the matrix element 



dxe^ 



(0|T(7r'^(x)4^(0))|0) (C14) 



where translational invariance has been invoked and the pion momentum is outgoing. Near the pion pole, it becomes 

m 

dxe^^-^ (0| r(7r'^(x)4'^(0) |0) -> 'J^^ ^ (C15) 
where the transition amplitude ■M'^p is given by 

tM%p,x''{p,P) - br% (P' - ^l)] \p2^^2 ■ (C16) 

The latter is a consequence of the isovector axial current being dominated by the pion-resonance contribution near 
the pole (4d| . Moreover, for leptonic decay, the transition amplitude is usually parametrised as 

zM'Xp = if.P^. (C17) 

Putting everything together leads to the expression 

\iT% _ ml)] = ihmlx'^ip, P). (C18) 

The evaluation of the pion decay constant is thus reduced to that of the pseudo-scalar contribution near the pion 
pole. Consider first a theory without any vector or axial four-quark couplings. The pseudo-sclalar contribution near 
the pion pole is given at leading Nc order by 

^% « |?^y|(1-GpJpp(P))<5''%(P) 

+ 2ml J dkf^ik+)fik_)TT [j5T''S,ik+h5T'S,{k_)] I (C19) 

where the second term is due to finite current quark mass. Using the definition of the pion wavefunction, the pion 
decay constant is extracted from Eq. (jC18|) and reads 

UmlS'"' = -TO« J dkTi[r{k,P)Sg{k+)i-f5T'Sg{k.)] 

+ ^{l^GpJppiP))d^'ls{P). (C20) 

where at the pole the second term is zero. Note also, that because the quark-pion coupling scales like 1/VlVc, the 
pion decay constant will have a \/Nc dependence. 

Introducing vector and axial couplings, the expression for the pion decay constant then becomes 



/^m^J-" = -ml J dkTT[nk,P)S,{k+)i^,T''S,ik^)] 

+ |LA.(P)|/5(P) + -^/v.(P)|5«^ (C21) 



Again, only the first term survives at P^ — m^. 

Contrary to the electromagnetic decay of the J /tp into a dilepton, the pion decay constant does not dependent 
on the path. This is due to the fact that the pion couples to the divergence of the isovector axial current, i.e., its 
longitudinal part, which is entirely determined by the axial Ward identity. 
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APPENDIX D: SCATTERING AMPLITUDES FOR THE NON-LOCAL NJL MODEL 

1. J/tp +7r ^ D + D 



Ml, 

Ml, 



•Pp 



where t = (p,r — Pd)^ and u = (p^r — Pd)^- 



2. J/ip + n^D + D* 



(Dl) 
(D2) 
(D3) 



M^ 

Mr. 

MZ 
MZ 
M7f 



i 

KDHD')^iXAP^-PD^)F;iiD^^{t), 
KDiD')i^)'D%{P.-PD)Fl';,',^Q,){u), 



(D4) 

(D5) 
(D6) 
(D7) 
(D8) 
(D9) 



where t = {p„ — Pd*)^ and u = (p^ — Pjj)^. 



J/ip + D* +D* 



MT 



MZ' 

M^/ 
MZ' 
Mff" 

Mr/ 



it)v^;{p.-PD')F^ 

1/ 



i,D*(D")y''h 



F^'L,^„{u)V^; {p.-PD')F^Zm^,i-'^) 



7Pa 

TrD*(D') 
ppvp_ 



Kd' (5, ) (P^ - PD' ) Ff/, ,n,M 



ipD'(Di 
Pup 



where t = (p,r —po*) and u = (p,r — Pd' ) • 



(DIO) 

(Dll) 

(D12) 
(D13) 
(D14) 
(D15) 
(D16) 
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4. J/il) + p^D + D 



i 

= F%^D')^u)V^P iPp - Pd) F^^n^-^), (D20) 
■^4? = F%j,^{s,t), (D21) 



where t={pp- pn) and u = (p^ - p^) . 

5. J/il} + D + D* 



Mt' = Y.Ff6\D)it)n(Pp-Pi^')F;'kn)ii)^ (D24) 

^^6^' = E<5(..)W^»^(P^-PS)^^;S^(5)H' (D25) 

Mf/ = FX^^^^{t)V"^^,{j>p-pn^)F%^^^^{t), (D29) 

= (Pp - P5) F;%^^^^{u) (D30) 



where t= {pp- po' ) and u= {pp- pq) . 
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6. J/tp + p^D* + D* 



i 

i 

MZ/' = F^-^t^^,)iu)vS;{p,-Pn^)Ff^:^^,^{u), (D34) 
^ef = F^B^oMt), (D35) 
^r/' = KDHD.^i)'^^piPp~PD')K5HD.^i)^ (D36) 

MZg'' = F^D'iD,){u)V^h (PP - PD')KD'iDr^''^^ ^^^T) 

MZh' = E^Sm^^S^ (D38) 

i 



where t={pp- pn* ) and u = {pp-p^,) . 
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